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ABSTRACT 

Using the high sensitivity of XMM-Newton, we have studied the X-ray emission 
of the two proto-typical late-stage merger remnants, NGC 3921 and NGC 7252. In 
the case of NGC 7252, this is complemented by archival Chandra data. We investigate 
the nature of the discrete X-ray point source populations and the hot diffuse gas 
components in these two galaxies, and compare them in the light of their different 
merger ages and histories. 

We detect 3 candidate ultra- luminous X-ray point sources in NGC 3921 and at 
least 6 in NGC 7252, for which we have high spatial resolution Chandra data. These 
have luminosities ranging from ~1.4xl0 39 — 10 40 erg s _1 (for Ho = 75 km s" 1 Mpc -1 ). 
We expect these ULXs to be high mass X-ray binaries, associated with the recent star 
formation in these two galaxies. 

Extended hot gas is observed in both galaxies. We have sufficient counts in the 
XMM-Newton data to fit two-component hot plasma models to their X-ray spectra 
and estimate the X-ray luminosities of the hot diffuse gas components to be 2.75xl0 40 
erg s" 1 and 2.09xl0 40 erg s" 1 in NGC 3921 and NGC 7252, respectively. These lumi- 
nosities are low compared with the luminosities observed in typical mature elliptical 
galaxies (L^ ~10 41 ~ 42 erg s _1 ), into which these merger remnants are expected to 
evolve. We do not see evidence that the X-ray halos of these galaxies are currently 
being regenerated to the masses and luminosities seen in typical elliptical galaxies. 
The mass of atomic gas available to fall back into the main bodies of these galaxies 
and shock-heat to X-ray temperatures is insufficient for this to be the sole halo re- 
generation mechanism. We conclude that halo regeneration is most likely a long-term 
(>10 Gyr) process, occurring predominantly via mass loss from evolving stars, in a 
sub-sonic outflow stage commencing ~ 2 Gyr after the merging event. 

Key words: galaxies: individual: NGC 3921 - galaxies: individual: NGC 7252 - 
galaxies: evolution - X-rays: galaxies - X-rays: ISM - X-rays: binaries 



1 INTRODUCTION 

The hypothesis that elliptical galaxies may be formed from 
the merging of two disk galaxies was first proposed by 
Toomre & Toomre (1972). It is now known to be an es- 
sential part of galaxy formation, but there are few merging 
systems close enough to be studied in depth. 

The merger-remnant systems NGC 3921 and NGC 7252 
are the proto-typical products of major disk-disk merging, at 
the end of the merger evolutionary sequence (Read & Pon- 
man 1998). Their nuclei have already coalesced, and their 
optical surface brightness profiles are close to the r 1//4 de 
Vaucouleurs law seen in typical elliptical galaxies (Schweizer 



1996; 1982). However, they are not yet fully relaxed, and 
still display the disturbed morphology of post-merger ob- 
jects, with tails and loops, which are the classic signature of 
equal- mass mergers (Figure 1). 

These two objects are nearby, and have been well- 
studied across a range of wavelengths. However it is only 
with the recent advent of XMM-Newton and Chandra that 
it has become possible to study their X-ray emission in any 
detail. Our XMM-Newton observations, and in the case of 
NGC 7252, archival Chandra data, in addition to the exist- 
ing multi- wavelength observations, therefore provide a new 
opportunity to study the X-ray emission - from both ex- 
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tended and point-like sources - arising in post-merger, proto- 
elliptical galaxies. 

NGC 3921 and NGC 7252 have different physical prop- 
erties (see Table 1), most likely as a result of their different 
post-merger ages, which are ~0.7 and ~1.0 Gyr respectively 
(Schweizer 1996; Hibbard & Mihos 1995) and differences 
in their parent galaxies. NGC 3921 is believed to be the 
product of a gas-poor— gas-rich, Sa— Sc (or SO— Sc) merger 
(Hibbard & van Gorkom, 1996), whereas NGC 7252 is the 
product of gas-rich— gas-rich, Sc— Sc merger (Hibbard et al. 
1994). As NGC 7252 is more evolved than NGC 3921, its 
optical emission is more relaxed and symmetrical. The gas- 
rich nature of its progenitor galaxies has induced a vehement 
starburst, with more luminous young globular clusters aris- 
ing from this starburst than are seen in the gas-poor— gas- 
rich merger remnant NGC 3921 (see Table 1). These galaxies 
offer us an ideal laboratory to investigate how the different 
merger histories and merger ages of these galaxies are re- 
flected in their X-ray emission. Hence, we investigate both 
the extended hot diffuse gas components and X-ray point 
source (XRP) populations in these galaxies. 

Ultra- luminous X-ray point sources (ULXs), with X-ray 
luminosities in excess of 10 39 erg s _1 , were detected in nor- 
mal and star-forming galaxies with the Einstein Observatory 
and ROSAT (Fabbiano 1995, Roberts & Warwick 2000), but 
it is only with the advent of XMM-Newton and Chandra that 
these studies could be extended outside our local environ- 
ment. The precise nature of these ULXs is still open to de- 
bate. King (2002) suggests that, in addition to the possible 
population of ULX supernovae remnants (SNR) , there are in 
fact two classes of ULX arising from super-Eddington mass 
inflow onto compact objects: thermal-timescale mass trans- 
fer in high-mass X-ray binaries (HMXBs), and long-lasting 
transient outbursts in low-mass X-ray binaries (LMXBs). 

The donor stars of HMXBs are expected to be O— B 
stars, with masses ~10 M©, at the end of their main se- 
quence lifetime. This sets the timescale for the evolution of 
HMXBs. The main sequence lifetimes of these donor stars 
are short (<300 Myr), so HMXBs are associated with re- 
cent star-formation. The ULXs observed in old stellar pop- 
ulations, e.g. elliptical galaxies, must therefore be of the 
LMXB type (e.g. Galactic micro-quasars), which is both 
longer-lived and takes longer to evolve. Hence, it is HMXBs 
that we expect to dominate in the merger-remnants NGC 
3921 and NGC 7252, both of which have substantial recent 
star- formation (Schweizer et al. 1996; Miller et al. 1997). 

An alternative theory is that ULXs are fuelled by ac- 
cretion onto intermediate-mass black holes (IMBH, 10 2 — 10 4 
M©), formed either from the collapse of massive popula- 
tion III stars or from merging or accretion in dense stellar 
regions (e.g. Colbert & Mushotsky 1999; Ebisuzaki et al. 
2001). This model is proposed to eliminate the necessity for 
super-Eddington accretion that HMXB models require. 

Massive X-ray halos (up to 10 10 M© of hot gas) are ev- 
ident in elliptical galaxies (e.g. Forman et al. 1985; Fab- 
biano et al. 1992). However, Read & Ponman (1998) studied 
the X-ray evolution of merging galaxies, with eight galaxies 
arranged in evolutionary order, and found that late-stage 
merger remnants such as NGC 3921 and NGC 7252 are 
under-luminous in X-rays when compared with typical el- 
lipticals. Log [(Lx/erg s _1 / (Ls/L©)] ~ 30.5 for a typical 
elliptical galaxy (from the Einstein observations of Fabbiano 



et al. (1992)), yet our estimates of the luminosities of the hot 
diffuse gas components give log [(Lx/erg s _1 / (Lb/L©)] ~ 
29.84 and 29.61 for NGC 3921 and NGC 7252 respectively, 
across the same energy range (0.2-3.5 keV) as the Fabbiano 
et al. results. If, as is expected, these two galaxies are to 
eventually evolve into normal ellipticals, their X-ray halos 
must be rebuilt. 

Three hypotheses have been proposed for the regener- 
ation of the hot halo: the infall of tidally stripped cold gas; 
mass loss from strong stellar winds and the return of hot 
gas expelled from the main body of the remnant during the 
nuclear merger stage. The sensitivity of XMM-Newton offers 
us the best opportunity to investigate these scenarios. 

We have ~25 ksec XMM-Newton European Photon 
Imaging Camera (EPIC) observations of NGC 3921 and 
NGC 7252, together with a 28.6 ksec archival Chandra ob- 
servation of NGC 7252, with which we can study both the 
ultra-luminous discrete source population and the hot dif- 
fuse gas in these two galaxies. We adopt Ho = 75 km s _1 
Mpc -1 throughout this work. 

The paper is structured as follows: §2 describes the ob- 
servations and data reduction; in §3 and 4, images and spec- 
tral fitting for NGC 3921 and NGC 7252 respectively are 
presented; in §5, the ultra-luminous discrete source popula- 
tion and the hot diffuse gas in the two galaxies are discussed 
and compared. §6 presents a summary of our results. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 XMM-Newton observations: NGC 3921 and 
NGC 7252 

The details of the XMM-Newton observations of NGC 3921 
and NGC 7252 are listed in Table 2. Only the EPIC data for 
these two galaxies are discussed in this paper; there are too 
few counts for meaningful analysis of the Reflection Grating 
Spectrometer data. The XMM-Newton Science Analysis Sys- 
tem (SAS), version 5.4.1, was used to process the observation 
data files into calibrated events lists, clean up background 
flaring (with flag = #XMMEA_EM for the MOS data, and 
#XMMEA_EP for the pn data) and filter the events. For the 
MOS1 and MOS2 data, single, double and quadruple pixel 
events are selected, and for the pn data single and double 
events. The 'good times' remaining for each camera after the 
cleaning process are listed in Table 2. For both galaxies, the 
data were very clean, with only a small fraction requiring 
removal because of background flaring. 

SAS and HEASARC's FTOOLS version 5.2 were used 
to create background-subtracted, exposure-corrected, all- 
EPIC mosaiced images in four energy bands: 0.2—0.5, 
0.5-2.0, 2.0-4.5 and 4.5-10.0 keV. The exposure maps 
were created using eexpmap, and the background maps with 
esplinemap, following the prescription given in the SAS 
data analysis threads. The mosaiced images were simply- 
smoothed with a Gaussian of FWHM 6.25" (compared with 
XMM-Newton spatial resolution ~7") to produce the X-ray 
contour maps presented in this paper. Source-searching was 
carried out for each of the three EPIC instruments in all 
four energy bands simultaneously, using the SAS source- 
searching packages. 

Two-dimensional surface-brightness fitting was carried 
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NGC 3921 



NGC 7252 



Parent galaxies 

Age / Gyr since peri-galacticon 
Log (L B I Lq) 
Log (L FIR I Lq) 
Seo / Sioo 

Optical isophotes of main body 
Central structure 

Nucleus 

Young globular clusters 
Distance / Mpc 
kpc / " 

Galactic line-of-sight N H / xlO 20 



Sa-Sc (or SO-Sc) 
gas-poor + gas-rich 

10.76 3 
10.18 4 
0.52 5 
'sloshing' 
young stellar populations are 
asymmetrically distributed 
LINER 
-9 < M v < -13 
77.9 
0.360 
1.14 6 



Sc-Sc 
gas-rich + gas-rich 
~1.0 2 
10. 81 4 
10.41 4 
0.57 5 
symmetrical 
equilibrated central disk of 
molecular and ionised gas 
quiescent 
-9 < M v < -16 
63.0 
0.304 



1.95 



Table 1. Summary comparison of the properties of NGC 3921 and NGC 7252 (for H = 75 km s -1 Mpc -1 ) l Schweizer 1996; 2 Hibbard 
& Mihos 1995; 3 Liu & Kennicutt (1995); 4 Fabbiano et al. (1992); 5 Read k, Ponman (1998); 6 Dickey k, Lockman (1990). 
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Figure 1. DSS optical images of NGC 3921 (left, R-band) and NGC 7252 (right, B^-band). The bright loop to the south and fainter 
tail to the northeast in NGC 3921, and the tails to the northwest and east in NGC 7252 are clearly visible. These loops and tails are 
typical of merger-remnant galaxies. 



out on the 0.3—1.2 keV mosaiced raw images, using the 
Sherpa software from the Chandra Interactive Analysis of 
Observations (CIAO) software package, version 2.3. Expo- 
sure correction is carried out on the model rather than 
the data, so as not to interfere with the fitting statistics. 
The background is modelled with a constant amplitude 
component, allowed as a free parameter in the fitting pro- 
cess. The area over which the fitting is carried out is small 
(37.4"x37.4" and 58.3"x40.7" for NGC 3921 and NGC 7252 
respectively, centred on the nucleus in both cases), so expo- 
sure and PSF effects on the background are negligible.. The 
image is convolved with the combined all-EPIC point spread 
function, weighted to account for the greater effective area 
of the pn camera compared with the MOS cameras. The 
best-fitting model is found using the Cash statistic. 

The hardness ratio maps of Figures 3 and 13 were 



calculated in the following way. First, the (unsmoothed) 
0.2-0.5, 0.5-2.0, 2.0-4.5 and 4.5-10.0 keV background- 
subtracted, exposure-corrected, all-EPIC mosaiced images 
described above had pixels with < lxlO -6 counts per sec- 
ond set to zero. The hardness ratios were then calculated 
using HR n = (E n+1 - E n )/(E n+1 + E n ) , where n = 1 - 3, 
Ei is the 0.2—0.5 keV image, E2 is the 0.5—2.0 keV image, 
E 3 is the 2.0—4.5 keV image and E 4 is the 4.5-10.0 keV im- 
age. The hardness ratio maps were then simply-smoothed, 
with a Gaussian of FWHM 6.25". Pixels with values ±1 
in these maps were then set to zero, to remove areas where 
photons are only detected in one band. For NGC 7252, there 
is no HR3 map, as there is so little emission in the hardest 
energy band. Within the outermost optical contours plotted 
in Figures 3 and 13, the uncertainties in the hardness ratios 



© 0000 RAS, MNRAS 000, 000-000 



4 L.A. Nolan et al. 




12 W D4 11S1:DD 12 IB D4 11:51 :DD 

Right aacenalon (hh:mmxfi.«) Right aflnenalon fhh:mm:«a«) 



Figure 2. DSS optical (R— band) images of NGC 3921 (greyscale with thin grey contours) with all-EPIC X-ray emission contours (thick 
black and white, as described in § 2) overlaid. Top left: 0.2—0.5 keV; top right: 0.5—2.0 keV; bottom left: 2.0—4.5 keV; bottom right: 
4.5—10.0 keV. The X-ray contour levels increase by factors of two from 1.375xl0 -6 counts s _1 . Crosses: globular clusters and stellar 
associations (Schweizer et al. 1996); triangles: NGC 3921 and the dwarf galaxy 2MASX J11505929+5504133 to the east (Two Micron 
All-Sky Survey Team, 2003. Final release). The peak of the nuclear X-ray emission is off-set from the centre of the optical emission. The 
dashed circle in the top right plot represents the half-energy width (14"). The pixel sizes are 2.5" by 2.5". 



are of the order a few percent for HRi and HR2, and ~10% 
for NGC 3921 HR 3 . 



2.2 Chandra observation: NGC 7252 

NGC 7252 was observed with the Chandra Observatory (ob- 
servation ID 2980) on the 17th November 2001 for approx- 
imately 28.6 ks, with the back-illuminated ACIS-S3 CCD 
chip at the focus. Data products, correcting for the motion 
of the spacecraft and applying instrument calibrations, were 
produced using the Standard Data Processing system at the 
Chandra X-ray Center (CXC). A preliminary analysis of 
these products was then performed using the CXC CIAO 
software suite (version 2.2.1). 

A light-curve extracted from a large non-source area 
revealed a moderate flare, with a sharp peak at ~twice the 
quiescent background level, occurring over the latter quarter 
of the observation. Therefore, we performed our analysis on 
the flare-subtracted data (which amounted to 22.15 ks). 



3 NGC 3921 
3.1 Spatial data 

The background-subtracted, exposure-corrected, smoothed 
X-ray contour maps of NGC 3921, prepared as described in 
§2.1, are presented in Figure 2, superposed on the Digitized 
Sky Survey (DSS) optical (R— band) image. Four energy 
bands are shown: 0.2-0.5, 0.5-2.0, 2.0-4.5 and 4.5-10.0 
keV. The strongest emission is seen in the 0.5—2.0 keV 
range, but even at high energies (4.5—10.0 keV), there is 
still significant emission, presumably, at least in part, from 
the low-ionization nuclear emission-line region (LINER) re- 
siding in this galaxy (Stauffer 1982a,b). The peak of the 
X-ray emission is offset to the south east of the optical peak 
(central triangle) , with the offset being greatest in the soft- 
est band. This is indicative of a possible absorption effect, 
obscuring the true X-ray nucleus in the soft band, which 
is consistent with the results of model fitting to the nuclear 
spectra discussed in §3.2. The X-ray emission is not extended 
beyond the main body of the galaxy. 

Figure 3 shows the three hardness ratio (HR) maps cal- 
culated for NGC 3921 as described in §2, with optical con- 
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object 


date 


camera 


obs time /ksec 


filter 


window 


good time / ksec 


NGC 3921 


27/04/2002 


MOS1 


26.0 


medium 


PrimcFullWindow 


24.8 






MOS2 


26.0 


medium 


PrimeFullWindow 


25.1 






pn 


23.7 


medium 


PrimeFullWindow 


20.5 


NGC 7252 


13-14/11/2001 


MOS1 


27.6 


medium 


PrimcFullWindow 


26.0 






MOS2 


27.6 


medium 


PrimcFullWindow 


26.2 






pn 


25.0 


medium 


PrimcFullWindow 


20.3 



Table 2. Details of XMM-Newton observations of NGC 3921 and NGC 7252. 



source 


RA(J2000) / h:m:s 


Dec(J2000) / d:m:s 


detection likelihood 


10 2 counts s 1 


3921-X1 


11:51:06.96 


55:04:42 


167.5 


8.913±0.272 


3921-X2 


11:51:07.92 


55:04:13 


56.7 


2.652±0.134 


3921-X3 


11:51:03.84 


55:04:50 


74.2 


0.365±0.081 


3921-X4 


11:51:07.68 


55:04:35 


56.5 


1.723±0.312 



Table 3. Sources detected in NGC 3921, using the SAS source-searching software. The positions of the sources are plotted in Figure 3. 



tours from the DSS (R— band) image superposed. In these 
plots, darker greys represent softer emission, and lighter 
shades are harder. The detected sources are marked with 
circles, and identified according to their listing in Table 3. 
Most emission is in the 0.5—2.0 keV range; the HRi map 
shows the hardness peaking (from ~0.35 in the outer regions 
to ~0.74) in the centre of the galaxy; for HR2 and HR3, the 
hardness ratio remains negative, but broadly declines in the 
inner regions by a factor ~2. 

These plots reveal an irregular, energy-dependent struc- 
ture to the X-ray emission. Although the central emission in 
the softest bands (HRi, top) is relatively symmetrical about 
the centre of the galaxy, at higher energies (HR2 and HR3), 
this is certainly not the case. The hardness peak in HR2 near 
the centre of the optical emission probably results from the 
LINER emission, which peaks in X-rays in the 2.0—4.5 keV 
band (see the spectrum in Figure 7). However, outside this 
very central region, the asymmetry may be due to irregular 
dust obscuration, unresolved X-ray point sources, or the dis- 
turbed distribution of the hot diffuse gas components. As is 
confirmed by the spectra plotted in Figures 7, 8 and 9, the 
LINER has much harder X-ray emission than the discrete 
sources, which have little emission above 5 keV. 

Interestingly, in Figure 2, we see what appears to be a 
'bridge' of X-ray emission, in the 0.5—2.0 keV band, from the 
main body of the remnant towards the dwarf galaxy 2MASX 
J11505929+5504133. Schweizer (1996) finds a recession ve- 
locity for this galaxy (S96 G4 in his notation) of 5712 km 
s~\ so close to that of NGC 3921 (5838 km s _1 ), that it is 
likely to be physically associated with it. The existence of 
the bridge is suggestive of the formation of the dwarf galaxy 
from tidal debris ejected during the merging of the parent 
galaxies, as may be the case at the end of the southern tail 
of the Antennae galaxies (e.g. Zwicky 1956; Deeg et al. 1998; 
Hibbard et al. 2001). However, the optical surface brightness 
of S96 G4 is much higher than that seen in any such candi- 
date dwarf galaxies (Schweizer 1996). The 0.5—2 keV emis- 
sion in the bridge region is ~ 4<r above the background emis- 



sion, and so this feature is likely to be real. In addition to the 
suggested 'bridge', there are two arms of Ha emission from 
the main body of NGC 3921 towards this galaxy (Hibbard 
& van Gorkom 1996), which correspond closely to groups of 
young stellar associations and candidate young clusters im- 
aged with the Hubble Space Telescope (HST; Schweizer et 
al. 1996). These are visible in Figure 4, above and below an 
imaginary line connecting the dwarf galaxy with the nucleus 
of NGC 3921 (the lower arm is right at the bottom edge of 
the frame). We tentatively suggest the possibility that these 
two 'ridges' of Ha emission, which are nearly orthogonal to 
the main tidal-tail structures, are wakes, produced by this 
small, dense disk galaxy having punched through some of 
the dense HI gas south of NGC 3921 at relatively high ve- 
locity (several hundred km s _1 ). 

Source-searching with the SAS software finds four 
sources in NGC 3921, although there are too few counts 
from any of the sources to fit the spatial extent. The mini- 
mum detection likelihood for inclusion of the source in the 
output list is set at 8, so that the probability that a ran- 
dom Poissonian fluctuation has caused the observed source 
counts is 3.35xl0~ 4 . 3921-X1, the most significant source, 
is coincident with the nucleus of NGC 3921 (RA(J2000) 
Ilh51m06.9s, Dec(J2000) +55d04m43s; Two Micron All- 
Sky Survey Team, 2003. Final release), where the LINER 
resides. There is also a bright source ~ 0.5' south of the 
nucleus (3921-X2), associated with the southern tail and a 
concentration of stellar associations and/or globular clus- 
ters. This source is visible in all the energy bands, and is 
softer than the central region (which contains the LINER). 
There is some harder emission to the south of the main 
peak of 3921-X2, observable in the HR3 map. 3921-X3, to 
the north east, is detectable as a slight bulge in the X-ray 
contours in Figure 2, and is visible in the hardness ration 
maps. 

Figure 4 shows the Ha+[NII] (KPNO 2.1m telescope) 
and HI (VLA) emission detected in NGC 3921 (Hibbard & 
van Gorkom 1996). There is a ridge of Ha emission, which 
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Figure 3. Hardness ratio maps for NGC 3921, with optical con- 
tours from the DSS R— band image overlaid. Top: HR i(soft); 
middle HR2 (medium) bottom: HR3 (hard). See §3.1 for details. 
Lighter grey represents harder emission. The cross locates the 
co-ordinates of NGC 3921 (RA(J2000) Ilh51m06.9s, Dcc(J2000) 
+55d04m43s; Two Micron All-Sky Survey Team, 2003. Final re- 
lease*), and the circles the X-ra.v sources detected in the p~a.la.xv 




Figure 4. Hot, warm and cold gas in NGC 3921. Light con- 
tours - all-EPIC X-ray emission, 0.5—2.0 kcV, as described in 
§ 2. White patches: Ho (Hibbard & van Gorkom 1996); grey 
scale: DSS R— band image; dark contours: HI (Hibbard & van 
Gorkom 1996); crosses - globular clusters and stellar associations 
(Schwcizcr ct al. 1996). 3921-X2 is located in the same region as 
the tails of Ha and HI. 



is associated with recent star formation, along the southern 
tidal tail. This ridge coincides with an extended ridge of HI 
in the same direction. 3921-X2 appears to be associated with 
this region of star formation and atomic gas in the southern 
tail, and clearly resides within this tidal tail. 

We model the surface brightness of the 0.3—1.2 keV, 
all-EPIC mosaiced image of the main body of NGC 3921, 
within a square of side 37.4" centred on the nucleus. The 
model is a /3-model, with an adjustable power-law, of the 
form: 

f(r)=A[l + (r/r () ) 2 ]- a , (1) 

where r is the core radius, A is the amplitude at r = 
0, and a = 3/3 — |. A constant background component is 
also included, and the model is convolved with the mosaiced 
point spread functions (PSF) of the three EPIC instruments. 
The energy is cut off at 1.2 keV in order to exclude the bulk 
of emission from the LINER, whilst extending sufficiently 
far beyond the energy range most affected by absorption 
(i.e. ~0.5 keV) to retain sufficient counts to constrain the 
fitted parameters. 

Table 4 lists the best-fitting parameters of this model 
for the main body of NGC 3921. Figure 5 presents the data, 
best-fitting model and residuals of the fit (data — model). 
The HI and Ha contours are superimposed on the plot of 
the residuals. There is an excess to the south east of the 
nucleus, which is not spatially coincident with any globular 
cluster candidate. This is plausibly emission from 3921-X4, 
which we know has soft (< 1.2 keV) emission (see §3.2). 
However, we have too few counts from this source to allow 
us to successfully disentangle it from the hot diffuse gas 
emission and investigate its spectrum. 

The power-law index, a is 2.93T0.08, and the core ra- 
dius is 4.43T0.09 kpc (H = 75 km s^Mpc -1 ). The pa- 
rameters of the spatial model should be regarded as simply 
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-0.6 -0.4 -0.2 0.2 0.4 0.6 




Figure 5. Spatial fitting to the 0.3-1.2 keV, MOS1, MOS2 and 
pn mosaiced and smoothed image (top) of the main body of NGC 
3921. Top: mosaiced image, with contour levels overlaid; mid- 
dle: best-fitting 2-D /3-model, plus constant background, with the 
same contour levels as for the image overlaid; bottom: residuals 
(data — model) plus colour bar, with units in counts. The Ha con- 
tours (black), HI contours (white) and candidate globular clusters 
and stellar associations (crosses; Schweizer et al. 1996) are super- 
imposed upon the residuals. The diagonal cross marks the optical 
centre of NGC 3921. Each plot covers exactly the same region, of 
size 37.4" x37. 4", centred on the nucleus. The model is only fitted 
to the soft (< 1.2 keV) emission in order to exclude contributions 
from the LINER whilst retaining sufficient emission (~0.5 keV) 
which is not strongly absorbed. See §3.1 for details of the fitting. 



object core radius / kpc a amplitude / cnts 

NGC 3921 4.43±0.09 2.93±0.08 5.42±0.20 



Table 4. The results of 2-D spatial fitting to the 0.3-1.2 keV, 
MOS1, MOS2 and pn mosaiced images of the main body of NGC 
3921. The fitted model is a /3-model, as described in §3.1. rrj is 
calculated for Ho = 75 km s — 1 Mpc~ 1 . Figure 5 presents the data, 
best-fitting models and residuals of the fits. 



describing the shape of the surface brightness distribution, 
rather than offering a physical interpretation of that shape. 

Figure 6 compares the modelled X-ray surface bright- 
ness distribution with the mean V-band radial profile. The 
hot diffuse gas emission has a very different profile from the 
optical profile, with a more flattened core, and a fairly uni- 
form decline from ~3 kpc. The raw X-ray data, together 
with the best-fitting model surface brightness distribution, 
are also shown in Figure 6. 



3.2 Spectral fitting 

For the extraction of the spectra, source and background 
regions were defined with the help of the DS9 display tool. 
The source regions are shown in Figure 7, and have radii 
19.6" (nucleus), and 12.0" (3921-X2 and 3921-X3). The 
background regions were defined within large circles on the 
same CCDs as the sources, in regions uncontaminated by de- 
tected sources. There is probably some contamination from 
nuclear emission in the off-nuclear sources, but, from Figure 
6, the surface brightness of the nuclear region has dropped 
to <10% of the peak value for the hot diffuse gas at the 
distances of 3921-X2 and X3, so the X-ray emission will be 
dominated by the off-nuclear sources. SAS tasks were used 
for the extraction, and to generate the response products. 
The spectra were grouped so that there were at least 20 
counts per bin. Model spectra were fitted via x 2 minimisa- 
tion using HEASARC's XSPEC software package, version 
11.2. All the models fitted include a component for line- 
of-sight Galactic absorption, which is frozen in the fitting 
process, as well as intrinsic hydrogen absorption, which is 
allowed to vary. It is the fitted intrinsic absorption which is 
quoted in Tables 5, 6, 10 and 11. The Galactic line-of-sight 
equivalent hydrogen column density is given in table 1. 



3.2.1 The nucleus 

The 0.5-10.0 keV MOS1, MOS2 and pn spectra of the cen- 
tral region of NGC 3921, totalling 1987 counts, are shown 
in Figure 7. The extraction region is shown in the same 
figure. The spectra were modelled simultaneously by a two- 
component MeKaL emission spectrum for hot diffuse gas 
(Mewe et al. 1985; Mewe et al. 1986; Kaastra 1992; Liedahl 
et al. 1995) plus an absorbed power-law component. In Ta- 
ble 5, the parameters of the best-fitting model are listed, 
together with a spectral plot of the various contributions 
from the model components. Two MeKaL components in 
the model results in a much better fit than a single compo- 
nent, as one might expect in this disturbed starburst system. 
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model McKaL 1 MeKaL 2 power-law total 

kT / keV 0.79 (0.73-0.85) 0.45 (0.37-0.55) 

abundance 0.19 (0.16-0.22) 0.03 (0.01-0.05) 

photon index 1.50 (1.45-1.55) 

intrinsic N H / 10 22 cm" 2 1.93 (1.73-2.14) 

log (L x / erg s- 1 ) 40.29 (40.27-40.37) 39.90 (39.65-40.05) 41.33 (41.32-41.37) 41.38 (41.36-41.43) 

X 2 / d.o.f 93.2 / 86 



Table 5. Results from simultaneous fitting to the MOS1, MOS2 and pn, 0.5—10.0 keV, spectra of the nucleus of NGC 3921 (sec Figure 
7). The model fitted has components for the Galactic extinction, two MeKaL components, and an absorbed power-law. Abundance 
is measured relative to solar. The luminosities quoted are unabsorbed, for Ho = 75 km s" 1 Mpc - 1 . The errors quoted are the 90% 
confidence limits on each single parameter of interest, and the hydrogen column density refers to the intrinsic absorption, which is fitted 
in addition to the fixed line-of-sight Galactic absorption. 




r / kpc 




Figure 6. Top: diamonds: radial profile (raw data, binned over 
circular annuli) of the diffuse gas in NGC 3921, plus error bars, 
over the range of the 2-dimensional fitting; solid line: model radial 
profile, using the best-fitting 2-dimcnsional model parameters, 
convolved with the PSF. Bottom: dashed line: model radial profile 
of the hot diffuse gas in NGC 3921, using the fitted parameters 
from the best-fitting 2-D spatial fitting (sec Table 4and §3.1). The 
dotted lines represent the uncertainties in the model fitting. Solid 
line: mean V-band surface brightness profile (Schweizer 1996). 
The surface brightness has been normalised to unity at r=0 in 
both cases. 



The MeKaL components (kT = 0.79 and 0.45 keV) con- 
tribute most of the energy between 0.5 and 2.0 keV. The 
harder end of the spectrum (~ 3 keV) is well-modelled by 
an absorbed power-law component, with a photon index of 
1.50, and intrinsic hydrogen column density, xih = 1.93xl0 22 
cm~ 2 . The dominant contributor to this luminous (log (Lx 
/ erg s^ 1 ) = 41.33) power-law component is probably the 
low-ionization nuclear emission region (LINER) residing in 
this galaxy. However, it is likely that there are also some un- 
resolved X-ray point sources which contribute to the hard 
end of the spectrum. The overall quality of the fit is accept- 
able. 



3.2.2 3921-X2: the southern source 

We fit three different model spectra in turn to the 0.5—10.0 
keV spectra of 3921-X2, the southern source in NGC 3921, 
extracted from the region shown in Figure 7. The MOS1, 
MOS2 and pn spectra are fit simultaneously, and have a total 
of 322 counts. The three models are i) an absorbed power- 
law, frequently used to approximate the X-ray emission from 
X-ray binaries (XRBs); ii) a MeKaL thermal plasma model, 
and iii) a disk blackbody model (e.g. Mitsuda et al. 1984; 
Makishima et al. 1986), which models the X-ray spectrum 
from a standard accretion disk, consisting of multiple black- 
body components. XMM-Newton does not have sufficient 
spatial resolution to determine whether this source is point- 
like or diffuse, so we fit all these models in an attempt to con- 
strain its nature spectrally. If this extremely bright source 
is point-like, then it must be a ULX, as it has Lx > 10 39 
erg s" 1 (King et al. 2002). 

The results from the spectral fitting are shown in Ta- 
ble 6. The EPIC spectra, with the best-fitting model spec- 
trum overlaid, are presented in Figure 8. The best-fitting 
model is the disk blackbody model, with the temperature at 
the inner disk radius, Winner, = 0.81 keV. This is consis- 
tent with what one would expect from a black hole binary 
in the high/soft state (Maccarone et al. 2003). Although 
fitting the disk blackbody model results in the best fit, it 
should be noted (from the minimum values of x 2 ), that nei- 
ther of the other two models can be rejected. However, the 
high temperature of the best-fitting MeKaL model argues 
strongly against a hot gas origin for this source; typical X- 
ray emitting diffuse gas temperatures in merging systems 
are ~0.2— 0.5 keV (Read & Ponman 1998), and even in viri- 



© 0000 RAS, MNRAS 000, 000-000 



XMM-Newton observations of the late-stage merger-remnant galaxies NGC 3921 and NGC 1252 9 



SOURCE 3921-X2 3921-X3 

counts 322 149 

POWER-LAW 

photon index 2.37 (2.20-2.58) 1.28 (0.90-1.46) 

intrinsic N H / 10 22 cm" 2 0.21 (0.16-0.26) 0.00 (<0.13) 

log (L x I erg s- 1 ) 40.35 (40.26-40.43) 40.02 (39.92-40.22) 

X 2 / d.o.f 13.7 / 13 

CASH statistic / no. of bins 120.6 / 133 
MEKAL 

kT / keV 2.18 (1.94-3.90) 63.3(13.3- >9.9) 

abundance 0.00 (<0.35) 0.00 (<546.3) 

log (L x I erg s" 1 ) 40.23 (40.20-40.46) 40.01 (38.99- >42.07) 

X 2 / d.o.f 12.4 / 12 

CASH statistic / no. of bins 120.5 / 133 
DISK BLACKBODY 

T in I kcV 0.81 (0.96-0.70) 2.50 (1.62-3.95) 

intrinsic N H / 10 22 cm~ 2 0.00 (<0.11) 0.00 (<0.05) 

log (L x / erg s" 1 ) 40.15 (39.89-40.44) 40.00 (39.28-41.68) 

X 2 / d.o.f 11.3 / 13 

CASH statistic / no. of bins 122.5 / 133 



Table 6. Results from simultaneous fitting to the MOS1, MOS2 and pn, 0.5-10.0 keV, spectra of 3921-X2 and 3921-X3 (see Figures 
8 and 9). The models fitted include a component for Galactic absorption. Abundance is measured relative to solar. The luminosities 
quoted are unabsorbed, for Ho = 75 km s~ 1 Mpc -1 , with limits calculated from the confidence limits on the best-fitting parameters. 
The errors quoted are the 90% confidence limits on each single parameter of interest, and the hydrogen column densities refer to the 
intrinsic absorption, which is fitted in addition to the fixed line-of-sight Galactic absorption. See §3.2 for discussion of the fitting. 



alised systems they only reach ~1.0 keV (Matsushita et al. 
1994; Rangarajan et al. 1995). 

3.2.3 3921-X3: the north-western source 

We fit the same three models that we fit to 3921-X2 to the 
3921-X3 EPIC spectra. 3921-X3 has fewer counts than 3921- 
X2 (149 compared with 322), so for this source, the spectra 
were not grouped, and the Cash statistic was used, rather 
than x 2 minimisation. The Cash statistic does not provide 
an estimation of goodness-of-fit, but we can still compare 
the relative goodness-of-fit of the different models. 

The source extraction region, observed spectra and 
best-fitting model spectra are presented in Figures 7 and 
9. The best-fitting parameters for each model are listed in 
Table 6. The best-fitting model is the power-law model, with 
no intrinsic absorption, and an index of 1.28. The other two 
models cannot be statistically ruled out. However, the ex- 
tremely high temperature of the best-fitting MeKaL model 
(63.3 keV) argues very strongly against a hot gas origin for 
this source. It has a luminosity in excess of 10 40 erg s" 1 , so, 
if this is indeed a point source, it must be a ULX. 



4 NGC 7252 
4.1 Spatial data 

In Figure 10, the background-subtracted, exposure- 
corrected, smoothed X-ray contours for NGC 7252 in four 
different energy bands (0.2-0.5, 0.5-2.0, 2.0-4.5 and 4.5- 
10.0 keV) are presented, overlaid on the optical (Bj— band) 
DSS image. Globular and star clusters are marked with 



crosses. The X-ray emission is fairly symmetrical, with no 
shift in emission with energy. The optical and X-ray nuclei 
are well-aligned. The globular / star clusters are centrally 
concentrated. There are more of them and they are more 
luminous than those in NGC 3921 (Miller et al. 1997 and 
Table 1). There is negligible emission above ~ 4.5 keV. The 
central hot diffuse emission is not very extended in any of 
the energy bands. 

The CIAO tool wavdetect was used to search for point- 
like sources in the 28 ksec Chandra archival observation. 
Scales from 1 — 16 pixels (0.5—8") were investigated. In ad- 
dition to the extended nuclear emission, a total of 6 bright, 
high- significance (> 4a) point (i.e. small-scale) sources were 
detected in the 0.3—7 keV band within or close to the opti- 
cal confines of NGC 7252. Seven lower-significance (between 
2.3cr and 4a) point sources were also detected. A further 
source (X2) is obvious close to the bright nucleus, though 
not formally detected by the software (probably because of 
its proximity to the nucleus). Properties of this source have 
been estimated by calculating by hand the counts within 
an ellipse around this source, with the size and shape of 
the ellipse estimated by comparison with the surrounding 
sources: 7252-X1, X3, X9 and X10. The wavdetect intensi- 
ties were checked by performing aperture photometry us- 
ing custom-defined circular apertures of 1—8 pixel radius 
(1 pixel ~ 0.5"). The wavdetected counts were, within the 
uncertainties, equal to those obtained within small (~ 1—2 
pixels) radii, except for the nucleus, which is extended and 
therefore required a larger (~ 4.5 pixels) radius. 

The point-like nature of the sources is confirmed by 
comparing the Chandra radial profiles with the PSF, as 
shown in Figure 11. Here, the radial profiles of i) the nu- 
cleus, ii) the mean of the sources 1 — 14 (except 7252-X8, 
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Figure 10. DSS optical (Bj— band) images of NGC 7252 (greyscale with thin grey contours) with all-EPIC X-ray emission contours 
(thick black and white, as described in § 2) overlaid. Top left: 0.2—0.5 keV; top right: 0.5—2.0 keV; bottom left: 2.0—4.5 keV; bottom 
right: 4.5—10.0 keV. The X-ray contour levels increase by factors of two from 1.125xl0 -6 counts s _1 . Crosses: globular clusters and star 
clusters (Miller et al. 1997); triangles: galaxies (Maddox et al. 1990). The peak of the nuclear X-ray emission is aligned with the centre 
of the optical emission. Three off-nuclear sources are clearly detected, to the north-west, south-south-west and south. There is very little 
X-ray emission detected above ~ 4.5 keV. The dashed circle in the top right plot represents the half-energy width (14"). The pixel sizes 
are 2.5" by 2.5". 



which is not detected by Chandra), and iii) a 3 keV on- 
axis ACIS PSF, created using CIAO-mkpsf and scaled to 
the mean off-nuclear source profile, to which the estimated 
background rate is added. We take the mean profile of the 
off-nuclear sources, as the individual count rates are too low 
to justify comparison with the PSF. The mean non-nuclear 
source profile follows the PSF distribution and is therefore 
point-like. The nucleus, however, is definitely extended, out 
to 5" and beyond. 

The sources are listed in Table 7 and plotted in Figure 
12. 7252-X2 is the only source which may be associated with 
an observed globular or star cluster. 7252-X5 is a foreground 
star, and we exclude it from further discussion. 

In the XMM-Newton data, using the SAS source- 
searching software, we identify, in addition to the nucleus, 



five off-nuclear X-ray sources within the body of the galaxy 
from these data (see Table 8). The whole area covered by 
the Chandra observation was searched. The minimum detec- 
tion likelihood for inclusion of a source in the output list is 
10, so that the probability that a detected source is random 
Poissonian noise is 4.5x10 . 

None of these sources is obviously associated with any 
observed globular or star clusters. Although the spatial 
resolution of XMM-Newton is not sufficient to determine 
whether these sources are point-like or extended, two of 
the five XMM-Newton sources (7252-X1 and X3) are con- 
firmed (with a > 4.0) as point sources by the Chandra de- 
tections. 7252-X7 is also identified as point source by the 
Chandra observation, but with lower significance (a = 2.7). 
The co-ordinates of these sources for which there are de- 
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channel energy (keV) 



Figure 8. Spectral fit to 3921-X2, the southern source in NGC 
3291. Crosses: pn (upper) and MOS (lower) spectra; solid lines: 
best-fitting disk blackbody model spectra. See Table 6 for details 
of the models and §3.2 for description of the fitting and discussion. 




Figure 7. Spectral fit to the nucleus of NGC 3291. Top: spectral 
plot of this best-fitting model (top line) showing the two MeKaL 
components, which dominate at lower energies, and the substan- 
tial absorbed power-law component, which dominates at energies 
~ 2 keV. Middle: crosses: pn (upper) and MOS (lower) spectra; 
solid lines: best-fitting model spectra. Bottom: X-ray contour map 
(all-EPIC mosaic, 0.5—2.0 keV) showing the central region from 
which the spectra were extracted. See Table 5 for details of the 
models, and §3.2 for description of the fitting and discussion. 




channel energy (keV) 



Figure 9. Spectral fit to 3921-X3, the north-western source in 
NGC 3291. Crosses: pn (upper) and MOS (lower) spectra; solid 
lines: best-fitting power-law spectra. See Table 6 for details of the 
models and §3.2 for description of the fitting and discussion. The 
spectra in the plot have been rebinned to at least 10 counts per 
bin for clarity. 



tections in both the XMM-Newton and Chandra data agree 
very well. However, XMM-Newton source 7252-X4+ lies be- 
tween 7252-X4 and 7252-X12. It most likely includes con- 
tributions from the (low-significance) sources 7252-X12 and 
7252-X11. Hence, we label it 7252-X4+ to differentiate it 
from the Chandra source 7252-X4. 

This brings us to a total of seven significant ultra- 
luminous off-nuclear X-ray sources in NGC 7252 (7252- 
XI— X8, excluding 7252-X5), six of which we know to be 
point-like. There are in addition a further six point sources 
detected only in the Chandra data with lower (2.3< a <4) 
significance. 

Figure 13 shows the two hardness ratio (HR) maps cal- 
culated for NGC 3291 as described in § 2, with the contours 
from the DSS optical (Bj— band) image superposed. In these 
plots, darker shades represent softer emission, and lighter 
shades are harder. It should be noted that the grey scales in 
each plot are not the same. There is no HR3 map because 
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source 


RA(J2000) / h:m:s 


Dec(J2000) / d:m:s 


significance 


counts 


log (L x / erg s v ) 


NUCLEUS 


22:20:44.76 


-24:40:41 


25.0 


256.89±18.5 


40.75 (40.71—40.78) 


7252-X1 


22:20:41.15 


-24:40:07 


8.8 


1 8 54+4 36 


39.60 (39.49—39.70) 


7252-X2 


22:20:44.23 


-24:40:43 


~6.7 


15.0 ±4.0 


39.52 (39.38—39.62) 


7252-X3 


22:20:44.32 


-24:41:04 


15.5 


34 27+5 92 


39 88 C39 79—39 941 


7252-X4 


22-20-44 36 


-24:41:47 


6.7 


15.15+4 


39 52 f39 51 39 531 


7252-X5 


22-20-45 70 

£4 X4 • X4\) • t J • t \ ) 


-24-40-05 


8.9 


91 fi7+4 «n 
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7252-X6 


22-20-48 36 

4U X4 • X4\) • lu • * > V 1 


-24:41:42 


5.0 




39 3fi f39 20 39 481 


7252-X7 


22:20:40.67 


-24:39:30 


2.7 


5.58+2.45 


39.08 (38.83-39.24) 


7252-X8 












7252-X9 


22:20:49.23 


-24:40:33 


2.3 


4.73+2.23 


39.00 (38.73-39.18) 


7252-X10 


22:20:46.03 


-24:41:13 


2.9 


5.74+2.45 


39.11 (38.86-39.25) 


7252-XU 


22:20:45.98 


-24:41:49 


2.4 


4.74+2.24 


39.00 (38.73-39.18) 


7252-X12 


22:20:44.52 


-24:41:56 


3.3 


7.33+2.83 


39.20 (38.99-39.34) 


7252-X13 


22:20:45.79 


-24:39:41 


3.2 


6.59+2.65 


39.15 (38.93-39.30) 


7252-X14 


22:20:51.34 


-24:39:11 


2.8 


5.72+2.45 


39.11 (38.85-39.25) 



Table 7. Sources detected in the Chandra NGC 7252 data set, using wavdetect. The positions of the sources are plotted in Figure 12. 
Sources detected with greater than 4<r significance are printed in bold type. The luminosities are calculated in the 0.5— 10.0 keV range, 
using a power-law model with an index of 1.5, and Galactic line-of-sight absorption. The unabsorbed luminosities arc quoted. Ho = 75 
km s -1 Mpc _1 . 7252-X5 is a foreground star. 



source 


RA(J2000) / h:m:s 


Dec(J2000) / d:m:s 


detection likelihood 


10 2 counts s 1 


NUCLEUS 


22:20:44.79 


-24:40:43 


5109.8 


6.392+0.173 


7252-X1 


22:20:41.20 


-24:40:07 


108.7 


0.499+0.071 


7252-X2 










7252-X3 


22:20:44.45 


-24:41:06 


144.1 


0.858+0.091 


7252-X4+ 


22:20:44.48 


-24:41:51 


83.9 


0.461+0.068 


7252-X5 










7252-X6 










7252-X7 


22:20:40.79 


-24:39:31 


26.2 


0.227+0.060 


7252-X8 


22:20:46.41 


-24:40:48 


15.9 


0.300+0.073 



Table 8. Sources detected in the XMM-Newton NGC 7252 data set, using the SAS source-searching software. The positions of the 
sources are plotted in Figure 13. Sec also Figure 12 for the CTianaira-idcntified sources. 7252-X5 is a foreground star. 



there are too few photons detected in the hardest energy 
band. Most emission is in 0.5—2.0 keV range; the HRi map 
shows the hardness peaking in the centre of the galaxy. In 
HR2 the hardness ratio declines towards the centre. The 
positions of the five (XMM-Newton) off-nuclear sources are 
over-plotted on the HRi map. 

Figure 14 shows the Hq+[NII] (KPNO 2.1m telescope) 
and HI (VLA) emission detected in NGC 7252 (Hibbard et 
al. 1994). Hq+[NII] is mostly concentrated in the nucleus, 
with a loop extending out to the west. HI is associated only 
with the tails, both to the north-west and to the east. None 
of the off-nuclear X-ray sources is associated with the Ha 
region, except 7252-X2 (Figure 12), which lies within the 
central Ha region. The lop-sidedness of the HI emission is 
not reflected in the X-ray emission. 7252-X1, X2 and X7 are 
spatially correlated with the region where the north-western 
HI arm comes across the front of the remnant body (Hibbard 
and van Gorkom 1996), but we cannot tell whether they lie 
within, in front of, or behind the HI region. 

We model the surface brightness distribution of the 
main body of NGC 7252, within a rectangle 58.3" by 40.7", 
centred on the nucleus, following the same prescription as 
for NGC 3291. The energy range is restricted to 0.3-1.2 



keV, in order to exclude contributions from the power-law 
component (see §4.2 and the spectrum in Figure 17). The 
central region hosts at least one point source (7252-X2), so 
this is masked out of the fit (see Figure 15). Table 9 lists the 
best-fitting parameters for the /3-model, and Figure 15 plots 
the data, best-fitting model and fit residuals (data — model) . 
Again, the HI and Ha emission contours are over-plotted 
on the residuals image. Unlike the marked correspondence 
in the main body of NGC 3921, in NGC 7252, there is no 
correlation between the spatial distribution of HI and the 
deficit or excess of X-ray emission. However, the X-ray ex- 
cess, masked out of the fit, is in the very central region 
where the Ha peaks. Candidate globular clusters and stellar 
associations (Miller et al. 1997) are marked on the residuals 
plot, but these too do not obviously correspond with ex- 
cess (or deficit) emission (except in the case of the masked- 
out central region). The power-law index, a = 1.31+0.02, 
is less steep than that for NGC 3921, and NGC 7252 has 
a much more compact core, with a core radius 1.05+0.03 
kpc. As with NGC 3921, this model should be interpreted 
as a description of the shape of the hot diffuse gas surface 
brightness distribution, rather than as a physical model of 
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NGC7252 - Chandra - Radial profiles 
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Figure 11. Radial profiles of i) the nucleus (solid line), ii) the 
mean of the sources 7252-X1— 14 (except 7252-X8, which is not 
detected by Chandra) dotted line, and iii) a 3 keV on-axis ACIS 
PSF, scaled to the mean off-nuclear source profile, to which the 
estimated background rate is added (dashed line). The mean off- 
nuclear source profile is point-like, whereas the nucleus is ex- 
tended beyond the PSF. 




Figure 12. Sources detected in the XMM-Newton (white circles) 
and Chandra (black ellipses) observations of NGC 7252, overlaid 
on the DSS (Bj) image. Large fonts refer to those sources detected 
with >4cr significance, and small fonts to those with 2.3< a <4 
significance. 7252-X5 is a foreground star. 



it. Source 7252-X8 is most likely responsible for the excess 
of emission to the east and slightly south of the nucleus. 

Figure 16 compares the modelled X-ray surface bright- 
ness distribution with the mean V-band radial profile. The 
hot diffuse gas emission has a compact central core, al- 
though, as for NGC 3921, it is not as compact as the central 
V-band radial profile. The X-ray surface brightness distribu- 
tion has a similar slope to that of the optical surface bright- 




Figure 13. Hardness ratio maps for the all-EPIC observation of 
NGC 7252, with optical contours from the DSS image (Bj— band) 
overlaid. Top: HRi (soft); bottom: HR2 (medium). See § 2 for 
details of the hardness ratio calculation. Lighter grey repre- 
sents harder emission. The cross locates the centre of NGC 7252 
(RA(J2000) 22h20m44.775s, Dcc(J2000) -24d40m41.83s; Miller et 
al. 1997), and the circles the X-ray sources detected in the galaxy, 
which are labelled in the top plot as in Table 

ness distribution at larger radii (~4 kpc). The raw X-ray 
data, together with the best-fitting model surface brightness 
distribution, are also shown in Figure 16. 

4.2 Spectral fitting 

4.2.1 Nucleus 

The 0.5-10.0 keV MOS1, MOS2 and pn spectra of NGC 
7252, totalling 1375 counts, extracted from the region shown 
in Figure 17, were modelled simultaneously by a two- 
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model McKaL 1 McKaL 2 power-law total 

kT / keV 0.72 (0.66-0.77) 0.36 (0.33-0.41) 

abundance 0.14 (0.12-0.16) 0.64 (0.54-0.72) 

photon index 1.72 (1.66-2.39) 

intrinsic N H / 10 22 cm" 2 1.71 (1.26-3.66) 

log (L x I erg s- 1 ) 40.12 (40.07-40.15) 39.89 (39.80-39.98) 40.28 (40.06-40.31) 40.60 (40.47-40.64) 

X 2 / d.o.f 58.4 / 46 



Table 10. Results from simultaneous fitting to the MOS1, MOS2 and pn, 0.5-10.0 keV, spectra of the nucleus of NGC 7252 (see Figure 
17). The model fitted has components for the Galactic extinction, two MeKaL components, and an absorbed power-law. Abundance 
is measured relative to solar. The luminosities quoted arc unabsorbed, for Ho = 75 km s -1 Mpc - 1 . The errors quoted are the 90% 
confidence limits on each single parameter of interest, and the hydrogen column density refers to the intrinsic absorption, which is fitted 
in addition to the fixed line-of-sight Galactic absorption. See §4.2 for description of the fitting and discussion. 




Figure 14. Hot, warm and cold gas in NGC 7252. Thick contours: 
all-EPIC, 0.5—2.0 keV emission, as described in § 2. black patches: 
Ha (Hibbard et al. 1994); thin contours: DSS optical (Bj — band) 
image; white-on-grey scale: HI (Hibbard et al. 1994); crosses - 
globular clusters and star clusters (Miller et al. 1997). 



there is negligible emission above 4.0 keV. The goodness-of- 
fit of the best-fitting model is acceptable. The hot diffuse gas 
components in NGC 7252 are slightly cooler than those in 
NGC 3921. The luminosity of the nucleus estimated from the 
Chandra data set does not agree well with the XMM-Newton 
results, but this is a result of the different models fitted; fit- 
ting a power-law with a slope of 1.5 to the XMM-Newton 
data results in the same luminosity as for the Chandra data, 
within the uncertainties. We are able to fit a much more 
complex model to the XMM-Newton spectra, as there are 
far more counts available than in the Chandra data. This 
gives us an acceptable goodness-of-fit, which is not the case 
with the simple power-law model. 

Awaki ct al. (2002) speculate that the hard (2.0-10 
keV) X-ray emission that they detect in their 43.5 ks ob- 
servation of this galaxy using the Solid-State Imaging Spec- 
trometer on ASCA (1998) indicates the existence of nuclear 
activity or an intermediate-mass black hole in the remnant 
body. However, there is no other evidence of nuclear activ- 
ity in this well-studied object. The ASCA spectrum is ex- 
tracted from a circular region, 6' in diameter, centred on the 
nucleus. This region includes all of the sources we detect, so 
it seems more likely that the hard X-ray component in the 
ASCA spectrum represents the emission from unresolved X- 
ray point sources. The total X-ray luminosity observed by 
Awaki et al., 8.1xl0 40 erg s _1 , is consistent, within the un- 
certainties, with the sum of the luminosities of all the sources 
we detect. 



object core radius / kpc a amplitude / cnts 

NGC 7252 1.05±0.03 1.31±0.02 15.2±0.8 



Table 9. The results of 2-D spatial fitting to the 0.3-10.0 keV, 
MOS1, MOS2 and pn mosaiccd images of the main body of NGC 
7252. The fitted model is a /3-model, as described in §3.1. ro is 
calculated for Ho = 75 km s -1 Mpc _1 . Figure 15 presents the 
data, best-fitting models and residuals of the fits. See §4.1 for 
discussion. 



component MeKaL model plus an absorbed power-law com- 
ponent. The best-fitting parameters of this model are listed 
in Table 10. Figure 17 shows a spectral plot of the various 
contributions from the model components. The power-law 
component is much less luminous than in NGC 3921, and 



4-2.2 Off-nuclear sources 

We fit absorbed power-law and disk blackbody models to 
four of the five off-nuclear sources detected in our XMM- 
Newton observation of NGC 7252. The source regions are 
shown in Figure 18. The background regions were defined 
within large circles on the same CCDs as the sources, in 
regions uncontaminated by detected sources. There is likely 
some contamination from nuclear emission in the source re- 
gion for 7252-X3, but, from Figure 16, the surface brightness 
of the nuclear region has dropped to ~5% of the peak value 
for the hot diffuse gas at this distance, so the X-ray emission 
will be dominated by the off-nuclear sources. 

As with 3921-X4, 7252-X8 has too few counts to suc- 
cessfully disentangle it from the diffuse X-ray background 
spectrum and investigate its spectrum. We do not fit MeKaL 
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SOURCE 7252-Xl 7252-X3 7252-X4+ 7252-X7 

counts 221 109 264 79 

POWER-LAW 

photon index 1.38 (1.10-1.63) 1.98 (1.58-2.29) 1.88 (1.42-2.14) 2.06 (1.20-3.32) 

intrinsic N H / 10 22 cm~ 2 0.00 (<0.08) 0.00 (<0.06) 0.20 (0.05-0.33) 0.00 (<0.19) 

log (L x / erg s- 1 ) 39.89 (39.76-40.05) 39.54 (39.45-39.71) 39.89 (39.80-40.10) 39.14 (38.54-39.95) 

CASH statistic / no. of bins 184.4 / 204 109.4 / 101 241.4 / 232 64.1 / 78 

DISK BLACKBODY 

T in I kcV 1.70 (1.21-2.67) 1.24 (0.70-1.59) 1.30 (1.04-2.32) 0.42 (0.24-2.90) 

intrinsic N H / 10 22 cm- 2 0.00 (<0.07) 0.00 (<0.02) 0.00 (<0.11) 0.00 (<0.22) 

log (L x I erg s- 1 ) 39.80 (39.18-40.54) 39.56 (38.52-39.99) 39.73 (39.33-40.71) 38.88 (37.46-42.04) 

CASH statistic / no. of bins 187.0 / 204 124.4 / 101 240.3 / 232 66.6 / 78 



Table 11. Results from spectral fitting to the MOS1, MOS2 and pn spectra of the off-nuclear sources in NGC 7252 across the 0.5 — 
10.0 keV energy range (see Figure 18). The models fitted include a component for Galactic absorption. The luminosities quoted are 
unabsorbed, for Ho = 75 km s -1 Mpc -1 , in the energy range 0.5—10.0 keV. The errors quoted are the 90% confidence limits on each 
single parameter of interest, and the hydrogen column density refers to the intrinsic absorption, which is fitted in addition to the fixed 
line-of-sight Galactic absorption. 



models to the four spectra, as all four of these sources are 
confirmed point sources (in the case of 7252-X4+, poten- 
tially the superposition of two or three point sources), and 
therefore do not have a hot diffuse gas origin. The models 
were fitted across the 0.5—10.0 keV energy range. The spec- 
tra are not grouped into bins, so the Cash statistic is used 
for the fitting, as for 3921-X3. The extraction regions, ob- 
served and best-fitting model spectra are plotted in Figure 
18. 

Three of these four sources (7252-Xl, X3 and X7) are 
best-modelled by power-law models, with indices 1.38, 1.98 
and 2.06 respectively, comparable within the uncertainties, 
and with no intrinsic absorption. The fourth source, 7252- 
X4+, is better-modelled by a disk blackbody model, with a 
small (Nh = 2.0xl0 20 cm ~ 2 ) amount of intrinsic absorp- 
tion, and Winner = 1-30 keV. Although this may be con- 
taminated with contributions from 7252-Xl 1 and 7252-X12, 
the Chandra data shows that these two sources are less than 
half as luminous as 7252- X4, so the spectrum should still 
be dominated by contributions from 7252-X4. Of course the 
possibility that the X-ray emission from XI 1 and/or X12 
varies considerably over the timescale of the two observa- 
tions cannot be rejected, and one or both of these sources 
may have been more more luminous in the XMM-Newton 
data, observed four days before the Chandra data. 

All these sources are very luminous, with Lx > 10 39 erg 
s 1 , and therefore they are all ULXs, if they lie within NGC 
7252. Indeed, they are some of the most luminous ULXs ever 
observed. In §5.2, we discuss the number of these sources 
that we expect to be background objects. 

There are too few counts to fit complex models to the 
Chandra data, so we estimate luminosities using a power-law 
model with an index of 1.5 and no intrinsic absorption. The 
Ch an dra- derived luminosity agrees, within the errors, with 
the XMM-Newton luminosity for the source 7252-X7 (which 
is not well-constrained by the spectral fitting to the EPIC 
spectra). 7252-Xl and 7252-X3 have only a minor disagree- 
ment between the two observations. One might expect differ- 
ences between the two estimates due to the different models 
used, but differences may also arise from intrinsic variation 



in the source. This is discussed further in §5.2. The lumi- 
nosity of 7252-X4, estimated from the Chandra observation, 
lies within the uncertainty of the luminosity estimated using 
the disk-blackbody model fitted to the XMM-Newton spec- 
trum (the best-fitting model). However, the possibility that 
the XMM-Newton spectrum of 7252-X4+ is contaminated 
with emission from the possible Chandra sources 7252-Xl 1 
and 7252-X12 cannot be ruled out; the sum of the X-ray 
luminosities of 7252-X4, Xll and X12, estimated from the 
Chandra data, also lies within the uncertainties in the lumi- 
nosity determined from the XMM-Newton data. 

If these low-significance sources are real, and lie within 
NGC 7252, they, like the more significant sources, are ULXs, 
with L x > fO 39 ergs" 1 . 



5 DISCUSSION 

5.1 Starburst and dynamical ages 

The dynamical age, i.e. the age since the last peri-galactic 
passage that induced the merger, of NGC 3921 is inferred 
from its morphology. The nuclei of the progenitor galaxies 
have already merged, but long tails remain, and this, to- 
gether with the off-centre nucleus and 'sloshing' isophotes 
imply that the merger age is 700±300 Myr, which is a few 
central crossing times (t cr ~ 5xl0 7 yr, Schweizer 1996). In 
the case of NGC 7252, N-body simulations suggest that ~ 
580h~ 2 Myr have elapsed since peri-galacticon, i.e. ~ 1 Gyr, 
for H = 75 km s^Mpc" 1 (Hibbard & Mihos 1995). Its sym- 
metrical isophotes and equilibrated central disk of molecular 
and ionised gas are consistent with it being a few 10 s yr fur- 
ther evolved than NGC 3291. 

The age of the starburst in NGC 3921 has been es- 
timated by Schweizer et al. (1996) at 250-750 Myr, for 
Z— 1.0—0.2 Z©, from colour-magnitude (V— I, V) diagrams 
of the candidate globular clusters and stellar associations, 
and the stellar population evolutionary synthesis models of 
Bruzual & Chariot (1993). The results suggest that globu- 
lar clusters have been forming for most of the time since the 
merger began, and possibly even prior to peri-galacticon. In 
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Figure 18. Spectral fits with best-fitting models as listed in Table 11 for the off-nuclear sources in NGC 7252. Centre: X-ray contour 
map (0.5—2.0 keV) showing the regions from which the spectra were extracted. Surrounding plots: crosses: pn (upper) and MOS (lower) 
spectra; solid lines: best-fitting model spectra. See Table 11 for details of the models and energy bands, and §4.2 for description of the 
fitting and discussion. The spectra in the plots have been rebinned to at least 10 counts per bin for clarity. 



NGC 7252, Miller et al. (1997), again using (V-I, V) colour- 
magnitude diagrams for globular clusters and the Bruzual & 
Chariot models, estimate ages of 650—750 Myr for 1.0—0.5 
Z@, together with a central < 10 Myr population, most 
likely of O stars and stellar associations. 



5.2 Ultra-luminous X-ray point sources 

In NGC 3921, excluding the LINER, we find two sources 
for which we have measured intrinsic X-ray luminosities > 
10 39 erg s _1 , and which, if they are indeed point-like, and 
reside within the galaxy, must therefore be ULXs. 3921-X4 
is most likely ultra-luminous as well; the source-searching 
software assigns it ~5 times as many counts as 3921-X3 
(Table 3), which has L x = 8.32xl0 39 erg s" 1 . However, it 
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Figure 15. Spatial fitting to the 0.3-1.2 kcV, MOS1, MOS2 and 
pn mosaiced image (top) of the main body of NGC 7252. Top: 
mosaiccd image, with contour levels overlaid; middle: best-fitting 
2-D /3-model, plus constant background, with the same contour 
levels overlaid; bottom: residuals (data — model) and colour bar, 
where lighter colours are positive, and darker colours are nega- 
tive. The units are counts. The image and residual plots have 
been simply-smoothed with a Gaussian of FWHM 2.75", after 
the fitting has been carried out. Overlaid are globular clusters 
and stellar associations (crosses; Miller et al. 1997), Ha contours 
(dark lines) and HI contours (light grey lines). The diagonal cross 
marks the optical centre of NGC 7525. See §4.1 for details of the 
fitting. The circle marks the nuclear region excluded from the fit. 
Each plot covers the same region, of size 58.3"x40.7", centred on 
the nucleus. 
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Figure 16. Top: diamonds: radial profile (raw data, binned over 
circular annuli) of the diffuse gas in NGC 7252, plus error bars, 
over the range of the 2-dimcnsional fitting; solid line: model ra- 
dial profile, using the best-fitting 2-dimensional model parame- 
ters, convolved with the PSF. Bottom: dashed line: model ra- 
dial profile of the hot diffuse gas in NGC 7252, using the fitted 
parameters from the best-fitting 2-D spatial fitting (see Table 
9 and §4.1). The dotted lines represent the uncertainties in the 
model fitting. Solid line: mean V-band surface brightness profile 
(Schwcizcr 1982, Miller ct al. 1997). The surface brightness has 
been normalised to unity at r=0 in both cases. 

is not possible to disentangle 3921-X4 from the diffuse gas 
X-ray emission to perform spectral fitting to measure its X- 
ray luminosity without more counts and / or better spatial 
resolution. 

In NGC 7252, the XMM-Newton observations reveal 
four off-nuclear sources with Lx > 10 39 erg s _1 , and a 
fifth, unidentified source, which is potentially a ULX as well. 
Source searching of the November 2001 Chandra data finds 
five such sources with a significance > 4a, three of which are 
coincident with XMM-Newton sources, together with a low- 
significance detection of 7252-X7, also detected by XMM- 
Newton. If we include the low-significance sources, the num- 
ber of ULXs detected in NGC 7252 could be as high as thir- 
teen. Using the logN— logS relation of Tozzi et al. (2001), we 
can estimate the expected number of background sources 
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Figure 17. Spectral fit to the nucleus of NGC 7252. Top: Spec- 
tral plot of the best-fitting model (top line) showing the two 
MeKaL components, which dominate the model, and the modest 
absorbed power-law component, which peaks at ~ 2 keV. Mid- 
dle: crosses: pn (upper) and MOS (lower) spectra; solid lines: 
best-fitting model spectra. See Table 10 for details of the models. 
Bottom: X-ray contour map (0.5—2.0 keV) showing the nuclear 
region from which the spectra were extracted. See §4.2 for de- 
scription of the fitting and discussion. 



not physically associated with the merger remnants. We ex- 
pect to find 0—1 background sources in NGC 3921 in the 
XMM-Newton data, and the same in NGC 7252, over the 
areas shown in Figures 2 and 10 respectively. In the Chan- 
dra data, over the area shown in Figure 12, ~4.5±2 back- 
ground sources are expected at or above the low detection 
limit (2.3a). 

The difference in the number of sources detected by 
XMM-Newton and Chandra in NGC 7252 may be due in part 
to time variability of the sources. The Chandra observation 
was carried out only four days after the XMM-Newton obser- 
vation, so if variability does account for the differences in the 
source detection and modelled luminosities, it would have to 
occur on a timescale of ~ a few days. Fabbiano et al. (2003) 
find spectral variability in all the non-nuclear ULXs in the 
Antennae, consistent with accreting compact XRB models, 
but this is on timescales of two months to two years. How- 
ever, the orbital periods of HMXBs can be ~1 — 10's of days, 
and periodic fluctuations, which can increase the (hard) X- 
ray luminosity by a factor ~ 10, corresponding to the peri- 
astron passage of the compact object, have been detected by 
Laycock et al. (2003), using the Burst and Transient Source 
Experiment. This is consistent with our apparent detection 
of variability over a timescale of 4 days. 

Alternatively, the greater spatial resolution of Chandra 
may allow detection of sources which XMM-Newton can not 
resolve; this is almost certainly the case with 7252-X2. The 
XMM-Newton detection limit (for a point source modelled 
by a power-law with an index of 1.5) is ~3.3xl0 39 erg s _1 in 
NGC 3921 and ~1.7xl0 39 erg s" 1 in the less-distant galaxy 
NGC 7252. However, the lower significance limit, o =2.3 
represents a detection threshold, for the same model, of 10 39 
erg s~ in the Chandra observation of NGC 7252; we are not 
able to detect the low significance Chandra sources in the 
XMM-Newton data. 

In a recent study, Swartz et al. (2003) have carried out 
a preliminary analysis of 84 near-by galaxies in the Chandra 
archive, spanning the range of Hubble morphological types, 
in a search for ULXs. They detected 120 ULXs with Lx > 
10 39 erg s _1 in these galaxies, only nine of which have Lx 
> 10 40 erg s _1 . These reside in four galaxies, all in their 
'spirals & irregulars' class, which includes merging galaxies 
and merger-remnants. No elliptical galaxies in this survey 
have Lx > 5xl0 39 erg s _1 . They find that both the number 
and mean X-ray luminosity of ULXs increase significantly in 
merging and interacting galaxies, and also correlate with far 
infrared (FIR) luminosity. The number of ULXs per galaxy 
residing in galaxies which are in the process of merging is 
~5, and the mean X-ray luminosity is ~7xl0 39 erg s _1 . Al- 
though the source detection thresholds in the heterogeneous 
Chandra study and our XMM-Newton and Chandra data 
are different, so we cannot make a statistical comparison, 
the numbers and mean luminosities of the ULX candidates 
we find in NGC 3921 and NGC 7252 do not conflict with the 
preliminary results of Swartz et al.. They do, however, point 
out that the trend towards brighter ULXs in galaxies with 
higher numbers of ULXs may be a statistical fluctuation, 
and is quite weak. 

NGC 3921 and NGC 7252, therefore, contain some 
ultra-luminous X-ray sources. Most of those in NGC 7252 
are confirmed as point sources by the Chandra observation 
(Figure 12, Table 7). The spatial resolution of XMM-Newton 
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is not sufficient to confirm the point-like nature of 3921- 
X2, but this is the most luminous X-ray source detected 
in these observations. We can compare our results with the 
Chandra ACIS-S observations of the merging system NGC 
4038/4039 (the Antennae galaxies, Zezas et al. 2002). Eight 
point sources with Lx > 10 39 erg s _1 are detected in this 
system (assuming a 5 keV bremsstrahlung model and Galac- 
tic line-of-sight absorption). The brightest of these sources 
has a luminosity of 8.9xl0 39 erg s" 1 (corrected to Ho = 
75 km s _1 Mpc _1 ), measured in the 0.1 — 10.0 keV range. 
The brightest off-nuclear sources in our observations are 
3921-X2 and 3921-X3, which have X-ray luminosities Lx 
= 1.38xl0 40 and 1.05xl0 40 erg s -1 respectively. So, even 
in the narrower range (0.5—10.0 keV) across which we do 
our spectral modelling, these two sources are brighter than 
the brightest Antennae source. However, Miller et al. (2003) 
fit more complex, disk blackbody plus absorbed power-law 
models to the XMM-Newton spectra of a sub-set of four of 
these eight bright Antennae sources, and determine X-ray 
luminosities of 1.0— 3.7xl0 40 erg s _1 for them, in the energy 
range 0.3-10.0 keV. 

The best-fitting power-law models for the eight Anten- 
nae sources have indices ranging from 1.21 to 1.94, with 
intrinsic N H = 5.0 - 2.6 x 10 20 cm' 2 (Zezas et al. 2002), 
which is comparable with what we find in NGC 3921 and 
NGC 7252 for the sources best-fit by power-law models. 

Chandra observations of the interacting galaxy system 
NGC 4676 (the Mice), find five ULXs with Lx ~5xl0 39 
(Read 2003). This system is similar to, though less evolved 
than, the Antennae galaxies. Although in this more distant 
galaxy (D = 88 Mpc) source confusion is likely even with 
Chandra's spatial resolution, this is again consistent with 
the picture of enhanced numbers of very luminous XRBs in 
merging systems. 

3921-X2 and 7252-X4+ are best-modelled by disk black- 
body spectra, with Ti nner = 0.81 and 1.30 keV respec- 
tively; these sources are plausibly black hole systems in the 
high/soft state. For these models, kT ~ 1 keV(M/ M©)~ 1/4 
(e.g. Makishima et al. 2000), so these temperatures are con- 
sistent with the presence of stellar-mass black holes in these 
objects (~ 23.3 Mq and 3.5 Mq, respectively). However, it 
should be noted, as pointed out in §4.2.2, that the XMM- 
Newton spectrum of 7252-X4+ could be contaminated with 
emission from the sources 7252-X11 and X12, although these 
sources are only detected with low significance (~ 3a). Note 
also that the other spectral models are not ruled out by the 
fitting statistics. 

Observations of early-type galaxies suggest that 
40—70% of their ULXs are associated with young globu- 
lar clusters (Angelini et al. 2001; Maccarone et al. 2003). 
Therefore, we have used existing HST optical observations 
of globular clusters (NGC 3921: Schweizer et al. 1996; NGC 
7252: Miller et al. 1997) to look for the optical counterparts 
to the X-ray sources in these two galaxies. 

If the XMM-Newton coordinates (absolute astrometry 
—4— 5") of NGC 3921 are accepted without any further cor- 
rection, then the position of the X-ray nucleus is within 2" 
of the HST-optical nucleus, and 3921-X2 corresponds to the 
stellar associations 127 and/or 128 from Schweizer et al. 
(1996). The coordinates of these associations are listed in 
Table 12. There is no coincidence between 3921-X3 and any 
optical source, and the nearest optical sources to 3921-X4 



are a candidate globular cluster (GC 68) and stellar associ- 
ation (A 126), both lying > 5" from 3921-X4. However, at 
the distance of NGC 3921, stellar associations have effective 
radii which are typically <J 1", and for globular clusters, it 
is even less, ~ 0.03". The astrometric accuracy of HST is 
-1-3" (e.g. Brammer et al. 2002). 

If we accept the Chandra coordinates (absolute astrom- 
etry <1") for NGC 7252 without any further corrections, 
the X-ray and HST-optical positions of the nucleus agree 
to ARA/ADec = -0.3"/+0.4", and there is good agreement 
between the positions of 7252-X1, X3 and X7 in the XMM- 
Newton and Chandra data sets (Figure 12). The XMM- 
Newton detection of 7252-X4+ agrees less well with the po- 
sition of the Chandra detection, but this is probably due 
to confusion between the sources 7252-X4, Xll and X12. 
The good agreement means we can be confident in com- 
paring the X-ray source positions with optical observations. 
7252-X2 corresponds closely (within 2") to both the stellar 
association MSW97-45 and the globular cluster MSW97-9 
(Miller et al. 1997). The other sources have no such opti- 
cal coincidences (except for 7252-X5 which is a foreground 
star). The absolute astrometry of Chandra is <1". 

Hence, for these two post-merger galaxies, the fraction 
of candidate ULXs associated with globular clusters / stellar 
associations is less than 20% (assuming all these sources are 
actually associated with the merger remnants). 

If XRBs are created in globular clusters (White et al. 
2002, Maccarone et al. 2003), then it is possible that they 
could be ejected from their parent clusters. 3-body inter- 
actions are likely to be the dominant mechanism for this 
ejection. These could almost always eject XRBs with Mbh 
< 10 M© from globular clusters, and virtually inevitably 
from young stellar clusters, which have lower escape veloc- 
ities (Kulkarni, Hut & McMillan 1993; Sigurdsson & Hern- 
quist 1993). From Quinlan (1996), the velocity of the kick 
applied to a binary in a 3-body interaction, Vkick, is: 

Vk^k/kms- 1 ~ 4O(m/lOM ) 3/2 (M/lOOM Q r 1 , 

where m is the mass of the companion star, and M the 
mass of the compact object. So, for the ~ 10 kms -1 required 
to escape from a young stellar cluster (Coleman Miller & 
Colbert 2001), then if m ~ 20 M©, M can be as high as 
10 3 M . However, for a dense globular cluster, where Vkick 
£30 kms -1 (Webbink 1985), then even for m = 20, M < 400 
M . 

It is also possible that XRBs may be 'kicked' out of 
their parent globular cluster during the collapse of the com- 
pact object by which they are powered. Even at very mod- 
est velocities (v~30-100 km s" 1 , Cordes & Chernoff 1998; 
Fryer & Kalogera 2001, respectively), they could travel sev- 
eral hundred parsecs over the timescales suggested above 
(0.3-9 kpc in 10-300 Myr at 30 km s _1 ). The distance of 
3921-X4 from its nearest optical neighbours is ~1.9 kpc, so it 
is plausible that it has been ejected from one of GC 68 or A 
126. 3921-X2 and 7252-X2 are the only sources which corre- 
spond to one (or two) stellar associations. The 'kick' velocity 
from supernova collapse is inversely proportional to the mass 
of the system (Cordes & Chernoff 1998). This suggests that 
ULXs residing at any significant distance outside the dense 
regions in which they are most likely formed must either be 
stellar-mass black holes, or must not have been ejected via 
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the kick from supernova collapse. In the supernova kick sce- 
nario, therefore, 3921-X2 and 7252-X2, which lie closest to 
their putative parent globular clusters, are the most likely 
candidates for intermediate-mass black holes (IMBHs), with 
masses 10 2 i M/M Q i 10 4 , which could reach ULX lumi- 
nosities with isotropical sub-Eddington emission. 

It is difficult, therefore, to distinguish between an 
isotropic IMBH or beamed HMXB origin for ULXs on basis 
of estimating mass from positional considerations. Rigorous 
constraints on mass are needed from observations of the ra- 
dial velocities of optical counterparts, or gravitational waves. 

IMBHs could form either via multiple 
mergers/ accretion, or from the collapse of massive, negli- 
gible metallicity, Population III stars. Both these scenarios 
require a high stellar density for there to be a reasonable 
probability for interactions to occur, either for the multi- 
ple merging / accretion which forms an IMBH in the first 
case, or for the capture of a massive stellar companion in 
the second. Their existence in ULXs was invoked in order 
to avoid the necessity for super-Eddington emission from 
ULXs. However, King (2003), using Chandra observations 
of the Cartwheel galaxy, suggests that, even for IMBHs ac- 
creting from massive stars, the inferred production rate and 
low accretion efficiency would require super-Eddington ac- 
cretion. He postulates that, unless IMBHs accrete from an 
unknown, non-stellar source, with very specific properties, 
most, if not all, ULXs associated with star-formation must 
be HMXBs. The case in favour of an IMHB origin for ULXs 
is reviewed in detail in Coleman Miller & Colbert (2003). 

It is possible that some of the detected X-ray sources 
are in fact groups of less-luminous X-ray sources. The spa- 
tial resolution of Chandra (~0.5") corresponds to ~150 pc at 
the distance of NGC 7252 (63.0 Mpc). Unbound stellar as- 
sociations have radii of order twice this (Blauw 1964), while 
globular clusters, which are gravitationally bound, have half 
light radii of ~10 pc (Schweizer et al. 1996), much less than 
the resolution at this distance. 7252-X2, the only off-nuclear 
X-ray source in the Chandra data with a likely optical coun- 
terpart, could therefore be the sum of X-ray emission from 
multiple objects residing in a globular cluster or stellar as- 
sociation. 

The contribution of X-ray emission from normal O stars 
can be readily discounted. A typical solar metallicity O-type 
star emits ~3.5xl0 34 ergs -1 in the 2—10 keV range (Helfand 
& Moran 2001), so it would require ~6xl0 4 O-type stars to 
result in the 2—10 keV luminosity seen in this object. This 
is ~3 times the total number of O-type stars in the entire 
Milky Way, and they would all have to reside in a region 
of radius ~150 pc. The lower spatial resolution of XMM- 
Newton corresponds to ~ 2kpc, but the observed luminosi- 
ties of the off-nuclear sources detected with XMM-Newton in 
NGC 3921 could only be reached with ~ a few xlO 5 O-type 
stars. 

A more promising possibility for source confusion is a 
superposition of bright X-ray binaries. We can estimate the 
number of 10 38 - 10 39 ergs 1 sources we would expect to 
find in NGC 3921 and NGC 7252 by comparison with the 
well-studied Antennae galaxies. Fabbiano, Zezas & Murray 
(2001) find 23 sources in this energy range in the Anten- 
nae. Normalising this number by the B— band luminosity, 
which is well-correlated with the total point source lumi- 
nosity (Colbert et al. 2003), and by comparison with the 



Object RA(J2000) / DEC(J2000) / ~ V V- I 
h:m:s d:m:s 

3921-X2 11:51:07.63 +55:04:16.53 

A 127 11:51:07.67 +55:04:17.8 24.9 +0.05+0.18 

A 128 11:51:07.72 +55:04:16.3 25.2 +0.28+0.21 



Table 12. Stellar associations corresponding to 3921-X2 
(Schweizer et al. 1996). 



X-ray luminosity function in the Antennae galaxies, we ex- 
pect ~26 and 29 such sources in NGC 3921 and NGC 7252 
respectively, with total aggregate luminosities of ~4.2 and 
4.9xl0 39 ergs" 1 . 

One could imagine a large group of XRB sources asso- 
ciated with an exceptional star cluster, formed in the ex- 
treme conditions of a galaxy merger. Of our ULX sources, 
only 3921-X2 and 7252-X2 have potential cluster counter- 
parts. 3921-X2 has more than twice the X-ray luminosity 
estimated above for all the 10 38 — 10 39 ergs -1 sources in 
NGC 3921, and even for the lower luminosity 7252-X2, the 
majority of the expected bright XRB sources in NGC 7252 
would have to reside in a single region 5150 pc across. If 
the proposed optical counterparts to 3921-X2 or 7252-X2 
were truly exceptional, for example, if they hosted massive 
star- formation, then the high density of 10 38 — 10 39 ergs -1 
sources might be conceivable, but this is not the case (Miller 
et al. 1997; Schweizer et al. 1996). 

It seems most likely, therefore, that the candidate ULXs 
we observe in NGC 3921 and NGC 7252 are mostly sin- 
gle ultra-luminous objects. They are likely to be HMXBs, 
which are associated with recent star-formation. We see 
more ULXs, with higher luminosities, than are observed in 
galaxies without recent star formation (e.g. Swartz et al. 
2003; Irwin et al. 2003). O— B stars at the end of their main 
sequence lifetime are the companions to the compact objects 
(neutron stars or black holes) in HMXBs, and take ~10— 300 
Myr to reach this point in their evolution. This timescale is 
entirely plausible, given on-going star-formation since peri- 
galacticon in both these galaxies. 

5.3 X-ray halo regeneration scenarios 

Read and Ponman (1998), in a study of the X-ray evolution 
of on-going mergers using ROSAT data, find that the late 
stages of mergers are under-luminous in X-rays compared 
with typical ellipticals. Although when the two nuclei coa- 
lesce, massive extensions of hot gas are observed (up to 10 10 
M Q , e.g. in the Antennae, Arp 220, NGC 2623), after this 
peak, the X-ray luminosity decreases. In fact, in post-merger 
remnants, Lx / L_g is ~ 5 times less than that of a typical 
elliptical (e.g. Fabbiano et al. 1992). If major mergers are in- 
deed the progenitors of normal elliptical galaxies, the X-ray 
halo must be in some way regenerated following the merging 
incident. 

In this section, we use our XMM-Newton observations 
to investigate possible mechanisms for the regeneration of X- 
ray haloes in NGC 3921 and NGC 7252. As we have outlined 
in the introduction, there are three hypotheses proposed for 
rebuilding the hot halo: the infall of tidally stripped cold 
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gas, mass loss from strong stellar winds, and the return of 
hot gas expelled earlier in the merging process. 

5. 3. 1 Infall of tidally stripped cold gas 

The HI gas in the two merger remnants is clearly associated 
with the debris from their gas-rich parent-galaxies, i.e. the 
tails and western loop in NGC 7252 and the southern tail 
in NGC 3921. HI gas is not, however, seen to the north 
of NGC 3921 (Mhi/L r < 0.02 M /L , Hibbard & van 
Gorkom 1996). 

In NGC 7252, the HI gas is not obviously associated 
with X-ray emission. In NGC 3921, the picture is somewhat 
different. The extremely luminous source to the south of the 
nucleus (3921-X2) is clearly associated with the HI tail, Ha 
emission and a large group of star clusters, suggesting that 
it is related to recent star formation. 

The spatial resolution of XMM-Newton is not sufficient 
to determine whether this source is point-like, and the pos- 
sibility that this source is not an X-ray binary but a local 
region of shock-heated diffuse gas cannot be dismissed on 
the basis of the spatial data. However, it is unlikely that 
the region is being shock-heated via winds from close-by 
stellar clusters; although these can reach temperatures > 2 
keV (the temperature of the best fitting MeKaL model to 
3921-X2), their X-ray luminosities are several orders of mag- 
nitude lower (< 10 36 ergs -1 ) than in 3921-X2 (e.g. Stevens 
& Hartwell 2003). 

The possibility that 3921-X2 is a shock region due to 
the infall of cold gas still remains. However, we can estimate 
the upper limit to the post-shock temperature, and compare 
this with the temperature of the best-fitting MeKaL model 
of this source. From Cavaliere et al. (1997): 

kT 2 = Imy + lm, (2) 
3 8 

where kTi is the pre-shock temperature, kT2 is the post- 
shock temperature, m p is the proton mass, and v is the pre- 
shock velocity of the infalling cold gas. The upper limit of 
the gas velocity is the velocity it would gain from falling into 
the galaxy from rest at infinity. We take the dynamical mass 
of the galaxy to be 1.06xlO n M (Hibbard & van Gorkom 
1996), and assume the pre-shock cold gas temperature is 
negligible. Putting these assumptions into equation 2, we 
estimate the upper limit of the post-shock temperature to 
be ~0.96 keV (~10 7 K). This is well below the lower limit of 
the best-fitting MeKaL model (1.94 keV), and suggests that 
shocked gas is not the origin of 3921-X2. 

However, it should be noted that freshly shocked gas 
may not be in either thermal or ionization equilibrium, 
and this can alter the observed properties of the post- 
shock gas. Takizawa (1999) discusses the thermal equilib- 
rium timescale, noting that the ion temperature is initially 
much higher than the electron temperature in the post- 
shock region (since most of the incoming kinetic energy is 
in the ions), and that the two temperatures equilibrate on 
an electron-ion interaction timescale. However, this effect 
would lead to a low observed electron bremsstrahlung tem- 
perature, whereas we see a best fitting MeKaL temperature 
higher than the upper limit to the post-shock temperature. 
Moreover, the equilibration timescale at the temperature 
and gas density in the region of 3921-X2 (~0.4 Myr) is much 



shorter than the post-shock timescale, inferred from the time 
elapsed since star formation (several hundred Myr). 

Calculations of the thermal and dynamical evolution of 
a fast adiabatically expanding gas (Breitschwerdt & Schmut- 
zler 1999) show that a delay in recombination can lead to 
departures from ionization equilibrium. A full set of dynam- 
ically and thermally self-consistent equations is needed to 
calculate the timescale for recombination, which is beyond 
the scope of this paper. However, Breitschwerdt & Schmut- 
zler suggest that the dynamical timescale for a hot plasma 
(T> 10 6 K) is much shorter than the intrinsic timescales 
(e.g. for recombination, ionization etc.). The ionization non- 
equilibrium gives rise to the appearance of emission lines 
representative of higher temperatures than the actual gas 
temperature, and the over-estimation of the temperature if 
the spectrum is fitted with collisional ionization equilibrium 
models. However, for 3921-X2, the upper limit to the post- 
shock temperature is ~0.96 keV, and it is difficult to see 
how a temperature as high as the best fitting MeKaL tem- 
perature, 2.18 keV, could be mimicked by this process. In 
addition, we do not see any emission lines in the spectra 
(Figure 8) which we would expect from this scenario. This 
lack of emission lines is consistent with the very low abun- 
dance of the best fitting MeKaL model. Hence, we conclude 
that the high temperature of the best fitting MeKaL model 
is inconsistent with a shock-heated gas model. 

In NGC 7252, we do not see an X-ray emission region 
analogous to that of 3921-X2, i.e. with associated clusters, 
HI and Ha emission. However, in Figure 14, we see a 'cav- 
ity' in the 0.5—2.0 keV X-ray emission to the north, which 
appears to correspond to the north-western HI tail where it 
crosses the hot gas. Since the tail comes across the front of 
the remnant body (Schweizer 1982; Hibbard & van Gorkom 
1996), the increased HI surface density in the cavity region 
(15-20 xlO 19 cm- 2 , figs. 5c and 5d in Hibbard et al. 1994) 
strongly suggests that the 'cavity' is due to X-ray absorp- 
tion by the neutral gas of the tidal tail in the foreground. 
This effect may of course be enhanced by the positions of 
7252-X1 and 7252-X7. 

If X-ray halos are indeed regenerated by infall along 
tidal tails, then maybe insufficient time since peri-galacticon 
has passed for substantial regeneration to have occurred (e.g. 
Hibbard & Mihos 1995). Considerable HI is evident out at 
least as far as the optical tails in NGC 7252 and beyond the 
optical loops to the south of NGC 3921. Hibbard & Mihos 
(1995) model the kinematics of the HI gas in NGC 7252 us- 
ing N-body simulations. They calculate that, of the HI origi- 
nally residing in the spiral progenitors of this object, at least 
3xl0 9 ft~ 2 Mq has already fallen inwards. However, the cen- 
tral regions of NGC 7252 are Hi-poor, so atomic gas falling 
into the remnant body from the tails must be converted to 
other forms, at a rate of 1—2 M Q yr _1 . The process of infall 
is a natural explanation for the conversion of HI to hot gas, 
stars and cold molecular gas. This hypothesis is supported 
by observations of both NGC 3921 and NGC 7252 which 
show the presence of all these three phases in the remnant 
bodies (e.g. Whitmore et al. 1993; Hibbard et al. 1994; Hib- 
bard & van Gorkum 1996; Schweizer 1996; Schweizer et al. 
1996; Yun & Hibbard 2001). 

Using the results of the Hibbard & Mihos simulations, 
we can estimate if there is sufficient HI still to fall back 
into the remnant body of NGC 7252 from the tidal tails, 
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to increase the hot gas mass (and hence luminosity) to that 
seen in a typical elliptical galaxy. 

The simulations indicate that ~6xl0 8 fo~ 2 M© of HI will 
return to within 5R e over the next 2>h~ 2 Gyr, i.e. one fifth 
as much again as has already fallen into the central regions. 
We estimate the mass of hot gas already present in the rem- 
nant body from the parameters of the MeKaL components 
fitted to the nuclear region, and assuming spherical symme- 
try. The mass is ~3xl0 8 Mq , within the spectral extraction 
region. 'Typical' X-ray bright elliptical galaxies (not includ- 
ing brightest group and massive cluster ellipticals, which 
may be considered a different class) have hot gas masses of 
1O 9 -1O 1O M (e.g. Forman et al. 1985). Even if all the HI 
gas still to return to the remnant body is converted to hot 
gas, with none involved in star formation or conversion to 
molecular gas, the total hot gas mass would be just below 
10 9 M Q , and it is unlikely that the conversion to hot gas is 
anywhere near this efficient. For comparison, ~3xl0 9 Mq of 
HI has already fallen back into the body of the galaxy, but 
we see only ~3xl0 8 Mq of hot gas, which, assuming all the 
hot gas has been converted from HI, with no contributions 
from stellar mass loss, represents an efficiency of only 10%. 
It is hard to see that HI infall in this galaxy would create 
sufficient hot gas to reach the mass seen in typical elliptical 
galaxies. 

We know that the northern tail of NGC 3291 is Hi-poor, 
but the southern tail contains ~ 2xl0 9 /i~ 2 Mq of HI, with 
the tail becoming more gas-rich along its length (Hibbard 
& van Gorkom 1996). However, there are no simulations of 
the kinematics of the atomic gas in NGC 3921 in the liter- 
ature, so we do not know what proportion of this gas will 
return to the remnant body. In addition, the irregular X-ray 
structure of NGC 3921 invalidates the assumption of spher- 
ical symmetry used to calculate the hot gas mass in NGC 
7252, so we cannot make the same calculation. However, as- 
suming its similar X-ray luminosity to NGC 7252 indicates 
a similar gas mass, this galaxy would require an additional 
mass of hot gas ~5xl0 8 Mq to take it to typical values. This 
only seems possible if the HI kinematics and/or conversion 
to hot gas in NGC 3921 are quite different from those in 
NGC 7252. 

Given the differences between the predicted hot gas 
masses of typical ellipticals and those predicted from our cal- 
culations involving HI masses, it is difficult to believe that 
HI infall can be solely responsible for the regeneration of 
the X-ray halo to luminosities typical of elliptical galaxies. 
Additional support for the existence of another mechanism 
for halo regeneration comes from the work of O'Sullivan et 
al. (2001a). They see a long-term trend of increasing X-ray 
luminosity with mean stellar age, extending over 15 Gyr, 
long after the HI tails seen in mergers have been completely 
assimilated into the body of the galaxy. If the X-ray lumi- 
nosity continues to rise after the completion of HI infall, 
then, even if HI conversion is responsible for the presence 
of some of the hot gas, there must be an additional source, 
which contributes to the hot gas mass after virialisation, for 
example, stellar mass loss, as discussed below. 

5.3.2 Stellar winds 

The models of Ciotti et al. (1991) investigate whether mass 
loss from stars could regenerate the X-ray halo. They pre- 



dict a super-sonic wind-driven flow in the early stages of 
post-merger elliptical evolution, arising from two sources: 
mass-loss from evolving stars, and secularly declining heat- 
ing from Type I supernovae. During this stage, a galaxy is 
X-ray faint, and the amount of hot X-ray gas, and hence the 
X-ray luminosity, decreases. This stage can persist for ~ 2 
Gyr, and is followed by a sub-sonic outflow. As the flow ve- 
locity decreases, this flow allows the build-up of hot gas, and 
results in rising X-ray luminosity over timescales ~10 Gyr. 
In this model, NGC 3921 and NGC 7252 would be in the 
first, wind-driven stage, and indeed, they are X-ray faint. 
The gas density and temperature are expected to fall off 
smoothly with radius, without the presence of shocks. We 
see no evidence that this is not the case (if we are correct 
in assuming that 3921-X2 is a ULX). However, in order to 
model the behaviour of gas density and temperature with ra- 
dius in detail, we would need both better spatial resolution 
than XMM-Newton can offer, and more counts. The sub- 
sonic outflow stage of this model, where X-ray luminosity 
increases gradually with time, is consistent with the obser- 
vations of elliptical galaxies discussed by O'Sullivan et al. 
(2001a). 

5.3.3 Infall from an outer reservoir of hot gas 

The third possible mechanism for the regeneration of the 
X-ray halo is the return of the hot gas which was expelled 
from the main body of the galaxy during the nuclear merging 
stage. This expelled gas can reach large radii before stalling. 
It then falls back into the galaxy, and shocks are expected 
as it collides with the existing halo gas. These outflows are 
typically asymmetrical (e.g. Read & Ponman 1998). We see 
little evidence for the regeneration of the X-ray halo via this 
mechanism in NGC 3921 and NGC 7252. We see little X-ray 
emission at large radii in either galaxy. This does not mean 
that we can reject the possibility that there is significant 
hot gas outside the main remnant bodies which is as yet too 
diffuse to have been detected in these observations. However, 
although the infall of hot gas may be delayed with respect 
to the cold gas, as the hot gas reaches larger radii than 
the cold, it will still occur on a galaxy dynamical time-scale 
(O'Sullivan et al. 2001b). As with the infall of cold gas, it is 
hard to reconcile this with the long-term (~10Gyr) trend of 
increasing X-ray luminosity with mean stellar age observed 
by O'Sullivan et al. (2001a; 2001b). 



6 CONCLUSIONS 

We summarise our conclusions as follows: 

The ULX populations 

In both NGC 3921 and NGC 7252, we observe ultra- 
luminous non-nuclear X-ray point sources. Where spectral 
fitting has been possible, the spectra are well-fitted by either 
power-law models (3921-X3, 7252-X1, 7252-X3, 7252-X7) 
or a disk blackbody accretion disk model (3921-X2, 7252- 
X4+). These energy distributions are typical of XRBs, and 
we suggest that these sources are HMXBs, associated with 
the recent star formation in these two galaxies. However, the 
possibility that the three off-nuclear sources in NGC 3921 
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are not point-like cannot be ruled out at present. The point- 
like nature of the four XMM-Newton sources in NGC 7252 
is confirmed by Chandra observations, which also reveal a 
further two XPSs, so that the total number of confirmed X- 
ray point sources observed in NGC 7252 (assuming all are 
associated with the remnant) is six. 7252- X8 has no corre- 
sponding Chandra detection, and hence cannot be confirmed 
as a point source. The differences in the evolutionary stages 
and progenitor galaxies of NGC 392f and NGC 7252 do 
not appear to be reflected in their ULX populations. The 
numbers and luminosities of the candidate ULXs observed 
in NGC 392f and NGC 7252 are compatible with the con- 
clusion of Swartz et al. (2003), that more luminous ULXs 
are found in interacting and merging galaxies than in other 
morphological types. 

The distribution of hot diffuse gas 

In NGC 7252, the more evolved merger remnant, the 
diffuse X-ray gas is more relaxed: it is centred on the opti- 
cal centre, is more symmetrical in the central regions and 
is cooler than the hot diffuse gas in NGC 3921. This is 
most likely a result of its greater dynamical age. Hot gas 
in NGC 7252 has a more compact core than that in NGC 
3921, reflecting differences between the progenitors of these 
two galaxies and their merger dynamics. When compared to 
NGC 3921, the more violent star formation in NGC 7252 is 
not reflected in a higher X-ray luminosity for the hot diffuse 
gas. 

X-ray halo regeneration 

Although substantial masses of HI are available to fall 
into the remnant bodies over the next several Gyr, it seems 
unlikely that this is sufficient to increase the mass of hot gas 
in either NGC 3921 or NGC 7252 to the masses (and hence 
luminosities) seen in typical ellipticals. We see no compelling 
evidence that massive X-ray halo regeneration is occurring, 
even in these late-stage merger remnants. It is likely that 
regeneration occurs over much longer timescales, of up to 
15 Gyr, predominantly via mass loss from evolving stars, 
in a sub-sonic outflow stage commencing ~ 2 Gyr after the 
merging event (Ciotti et al. 1991). 
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